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PHASELESS IMAGING BY REVERSE TIME MIGRATION: 
ACOUSTIC WAVES * 

ZHIMING CHENt AND GUANGHUI HUANG* 

Abstract. We propose a reliable direct imaging method based on the reverse time migration 
for finding extended obstacles with phaseless total field data. We prove that the imaging resolution 
of the method is essentially the same as the imaging results using the scattering data with full phase 
information. The imaginary part of the cross-correlation imaging functional always peaks on the 
boundary of the obstacle. Numerical experiments are included to illustrate the powerful imaging 
quality. 


1. Introduction. We consider in this paper inverse scattering problems with 
phaseless data which aim to find the support of unknown obstacles embedded in a 
known background medium from the knowledge of the amplitude of the total field 
measured on a given surface far away from the obstacles. Let the sound soft obstacle 
occupy a bounded Lipschitz domain D C M 2 with v the unit outer normal to its 
boundary T &. Let u l be the incident wave and the total field is u = u 1 + u s with u s 
being the solution of the following acoustic scattering problem: 


A u s + k 2 u s = 0 in 
u s = —u l on fp, 


2 \d, 


w-™) 


0 as 


+oo, 


( 1 . 1 ) 

( 1 . 2 ) 

(1.3) 


where k > 0 is the wave number. The condition G3) is the outgoing Sommerfeld 
radiation condition which guarentees the uniqueness of the solution. In this paper, 
by the radiation or scattering solution we always mean the solution satisfies the Som¬ 
merfeld radiation condition (jl.3D . For the sake of the simplicity, we mainly consider 
the imaging of sound soft obstacles in this paper. Our algorithm does not require 
any a priori information of the physical properties of the obstacles such as penetrable 
or non-penetrable, and for non-penetrable obstacles, the type of boundary conditions 
on the boundary of the obstacle. The extension of our theoretical results for imaging 
other types of obstacles will be briefly considered in section 4. 

In the diffractive optics imaging and radar imaging systems, it is much easier to 
measure the intensity of the total field than the phase information of the field mm 
nu. It is thus very desirable to develop reliable numerical methods for reconstructing 
obstacles with only phaseless data, that is, the amplitude information of the total 
field \u\. In recent years, there have been considerable efforts in the literature to solve 
the inverse scattering problems with phaseless data. One approach is to image the 
object with the phaseless data directly in the inversion algorithm, see e.g. mm- 
The other approach is first to apply the phase retrieval algorithm to extract the phase 
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information of the scattering field from the measurement of the intensity and then use 
the retrieved full field data in the classical imaging algorithms, see e.g. m- We also 
refer to PQ for the continuation method and pa Ham] for inverse scattering problems 
with the data of the amplitude of the far field pattern. In m some uniqueness results 
for phaseless inverse scattering problems have been obtained. 

The reverse time migration (RTM) method, which consists of back-propagating 
the complex conjugated scattering field into the background medium and computing 
the cross-correlation between the incident wave field and the backpropagated field to 
output the final imaging profile, is nowadays a standard imaging technique widely 
used in seismic imaging 0. In 012IS], the RTM method for reconstructing extended 
targets using acoustic, electromagnetic and elastic waves at a fixed frequency is pro¬ 
posed and studied. The resolution analysis in MIS] is achieved without using the 
small inclusion or geometrical optics assumption previously made in the literature. 

In this paper we propose a direct imaging algorithm based on reverse time mi¬ 
gration for imaging obstacles with only intensity measurement \u\ with point source 
excitations. We prove that the imaging resolution of the new algorithm is essentially 
the same as the imaging results using the scattering data with the full phase infor¬ 
mation, that is, our imaging function always peaks on the boundary of the obstacles. 
To the best knowledge of the authors, our method seems to be the first attempt in 
applying non-iterative method for reconstructing obstacles with phaseless data except 
[9] in which a direct method is considered for imaging a penetrable obstacle under 
Born approximation using plane wave incidences. We will extend the RTM method 
studied in this paper for electromagnetic probe waves in a future paper. 

The rest of this paper is organized as follows. In section 2 we introduce our RTM 
algorithm for imaging the obstacle with phaseless data. In section 3 we consider the 
resolution of our algorithm for imaging sound soft obstacles. In section 4 we extend our 
theoretical results to non-penetrable obstacles with the impedance boundary condition 
and penetrable obstacles. In section 5 we report several numerical experiments to show 
the competitive performance of our phaseless RTM algorithm. 

2. Reverse time migration method. In this section we introduce the RTM 
method for inverse scattering problems with phaseless data. Assume that there are 
N s emitters and N r receivers uniformly distributed respectively on T s = dB s and 
T r = where B s ,B r are the disks of radius R s ,R r respectively. We denote by Cl 
the sampling domain in which the obstacle is sought. We assume the obstacle D C £2 
and Cl is inside B s ,B r . 

Let rd(x,x s ) = 4>(x,x s ), where <fr(x,x s ) = ^H^\k\x — x s |) is the fundamental 
solution of the Helmholtz equation with the source at x s G T s , be the incident wave 
and \u(x ri x s )\ = \u s (x r ,x s ) + u l (x r ,x s )\ be the phaseless data received at x r G T r , 
where u s (x,x s ) is the solution to the problem nm-ijua with u l (x,x s ) = $>(x,x s ). 
We additionally assume that x s ^ x r for all s = 1,2, ...,7V s ,r = 1,2, ...,7V r , to avoid 
the singularity of the incident field u l (x,x s ) at r = x r . This assumption can be 
easily satisfied in practical applications. In the following, without loss of generality, 
we assume R r = rR s , r > 1. 

Our RTM algorithm consists of back-propagating the corrected data: 


A( x r , x s ) 


\u(x r ,x s )\ 2 - |tt 8 (a; r .,a; 6 )| 2 
u l (x r ,x 3 ) 


( 2 . 1 ) 


into the domain using the fundamental solution <fr(x r ,z) and then computing the 
imaginary part of the cross-correlation between u l (z,x s ) and the back-propagated 
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field. 


Algorithm 2.1. (RTM for Phaseless data) Given the data \u(x r ,x s )\ = 

| u s (x r ,x s ) + u l (x r ,x s )| which is the measurement of the total field at x r E T r when 
the point source is emitted at x s E T S; s = 1,..., 7V S; r = 1,..., N r . 

1° Back-propagation: For s = 1 7 ..., N s , compute the back-propagation field 

2ttR Nr 

v b (z,x s ) = -— -^2$(x r ,z)A(x r ,x s ), Vzea (2.2) 

r r =1 

2° Cross-correlation: For z £ Q, compute 

I( Z ) = —fc 2 Im . (2.3) 

It is easy to see that 

f fn \ 2 d 70 2 

/(z;) =-fc 2 Im |— — ^ r y^y]^(2:,a; g )j>(a; r ,2:)A(a; r ,a; 8 )|, Vzefi. (2.4) 

This is the formula used in our numerical experiments in section 5. By letting 
7V S , N r -A oo, we know that (|2.4|) can be viewed as an approximation of the following 
continuous integral: 


I(z) = — k 2 lm j J x s )$(x r: z)A(x r , x s )ds(x r )ds(x s ), V 2 : E 12. (2.5) 

We remark that the above RTM imaging algorithm is the same as the RTM 
method in [3] except that the input data now is A(x r ,x s ) instead of u s (x r , x s ). Hence, 
the code of the RTM algorithm for imaging the obstacle with phaseless data requires 
only one line change from the code of the RTM method for imaging the obstacle with 
full phase information. 


3. The resolution analysis. In this section we study the resolution of the 
Algorithm 2.1. We first introduce some notation. For any bounded domain U C M 2 
with Lipschitz boundary T, let ||r||^i(c/) = (||+ d^} 2 1| <f \\ 2 L2 ^) 1//2 be the 
weighted H 1 (U) norm and \\v\\ H i/ 2 ^ = (d^ 1 ||v||| 2 ( r ) + Ml p) 1 / 2 be the weighted 
j H 1 / 2 (T) norm, where djj is the diameter of U and 




\v{x) -v(y)\ : 


■ds(x)ds(y) 


1/2 


/r \x ~ y I 2 

By scaling argument and trace theorem we know that there exists a constant C > 0 
independent of dp such that for any (j) G C 1 (5), 


H<AIEv2 ( r D ) + W^/dv\\ H -xf 2 {TD) < Cmax(|0(a;)| + d D \V<t>{x)\). (3.1) 

xED 

The following stability estimate of the forward acoustic scattering problem is well- 
known m so]. 

Lemma 3.1. Let g E H rl / 2 (T^), then the scattering problem: 

Aw + k 2 w = 0 in M 2 \D, w = g on T^, 
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as r —y oo 


(3.2) 

(3.3) 







admits a unique solution w G H^ oc (M, 2 \D). Moreover, there exists a constant C > 0 
such that \\dw/du\\ H -i/ 2 {FD) < C\\g\\ H i/ 2 (rD) . 

The far field pattern w°°(x), where x = x/\x\ G -S 1 = {x G M 2 : \x\ = 1}, of the 
solution of the scattering problem (|3.2p - (|3.3|) is defined as (cf. e.g. j 8 j P. 67]): 


w°°(x) = 


i- r 
e *4 r 

V8nk Jt d 


der^ _ dwM ik& . y 
W[V) du(y) dv(y) 


ds(y). 


(3-4) 


It is well-known that for the scattering solution of (|3.2j) - (j3.3j) (cf. e.g. j3j Lemma 
3.3]) 


Im [ w^-ds = k [ \w°°(x)\ 2 dx. 

Jr D dv Js i 


(3.5) 


Now we turn to the analysis of the imaging function I(z) in (|2.5j) . We first observe 

that 


A(X r , X s ) 


u s (x r ,x s ) + 


\u ( x r , x s ) | 
u^ ( x r , ) 


u s {x r , x s )u L (x r ,x s ) 
u^ {x r , ) 


This yields 


(3.6) 


I(z ) = —/c 2 Im 


—/c 2 Im 
—/c 2 Im 



^>(z,x s )^>(x r , 

<L(2(,X s )<L(x r , 

<f>(z,x s )<f>(x r , 


z)u s (x r ,x s )ds(x r )ds(x s ) 

\u S (x r ,X s )\ 2 

z) — 77 - —ds{x r )ds{x s ) 

u l [x r ,x s ) 


?T(x r ,x s )?d(x r ,x s ) 

2 :) - 7 ---- ds(x r )ds (x s ). 

J ^(x r ,x s ) v v 


(3.7) 


The first term is the RTM imaging function with full phase information in [3| and 
thus can be analyzed by the argument there. Our goal now is to show the last two 
terms at the right hand side of dS3) are small. We start with the following lemma. 

Lemma 3.2. We have \u s (x r ,x s )\ < C(1 + A)d j D) 2 (/ci7 r ) _1 / 2 (/ci7 s ) _1 / 2 for any 
X'p G T^, g G Tg. 

Proof We first recall the following estimates for Hankel functions 0 ( 1 - 22 )- 
(1.23)]: 


, I». (,) «>I S (s) 


1/2 


7 it 


\/t > 0 . 


By the integral representation formula, we have 


s (x r ,x s ) = J ^ 


d®(x r ,y) du s (y,x s ) . . , . s 

u s [y,x s ) a ; , A - q - - ®(x r ,y) ) ds(y). 


dv(y) 


dv(y) 


(3.8) 


(3.9) 


By (|3.ip we have 

II*)ll J H' 1 /2(r £) ) + *)/^ zy ll//- 1 /2(r 0 ) < C(1 +■ kdD){kR r ) 1 ^ 2 . 

The lemma follows now from Lemma o and the fact that u l {y,x s ) = —$>(y,x s ) for 

y £ r D . □ 

Lemma 3.3. We have IH^^I > [2/(57re)]| In1 1 for any t G (0,1). 
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Proof. We use the following integral formula 

a —rt 


h, 


W(t) = r _—_ dr Vt > 0 

[) 7r J 0 r x / 2 (r-2i) x / 2 ’ ’ 


0 w ^ < -1/2 ( r - 2i)V2 l 

where Re (r — 2i) l/ ' 2 > o for r > 0. By the change of variable 


I H, 


(i) 


2 r° 

(01 > / 

^ J 0 


5 1 /2( s _2it)V2 


(is = 


2 r 
K Jo 


> 


s 1 / 2 \s — 2it\ 


[ , 6 

Jo 


\s — 2it\ + s 


ds 


7T j 0 \s — 2lt\ 
2 f 1 e~ s , 

> — / -ds, 

5tt J t s 


where in the last inequality we have used |s — 2it| < 5s for s > t. This completes the 
proof by noticing that f t s~ 1 e~ s ds > e _1 f t s~ 1 ds = e~ l \ lnf|. □ 

The following lemma gives an estimate of the second term at the right hand side 

of (E3- 

Lemma 3.4. We have 


k z I I <fr(z,x s )<fr(x r , z) ds(x r )ds(x s ) 

r s Jr r u l [x r ,x s ) 


< C(l + kd D ) 4 (kR s )~ 1/2 . 


Proof Denote = {(x r ,x s ) G Y r x T s : \x r — x s \ < l/(2fc)}. Let x r = 
R r (cos 6 r , sin 0 r ), x s = R s (cos6 s , sm6 s ), where 0 r: 0 s G [0, 27t] . Since R r = tR s , 
r > 1, we obtain easily 

\q _ Q | 

\x r — x s \ = R s y/l + r 2 — 2 r cos 10 r — | > 2i? s Vorsin r —— (3.10) 

by Cauchy-Schwarz inequality. Now for (x r ^x s ) G we have then either < #o 

or 7r — 6>o < 2 6>s ^ < 7r, where #o = arcsin yy G (0,7r/2). 

By (|3.8[) and Lemma 13.21 

[[ ®(z, x s )$(x r , z) — f r ’ ds(x r )ds(x s ) 

JJ n k u l {x r ,x s ) 

<C(l + kd D ) 4 (kR s )- 3/2 (kR r )~ 3/2 [[ \ ln(k\x r - x s \)\ds(x r )ds(x s ). (3.11) 

J J Ofc 

By (13.1011 we have 


[[ IKfcl 

J J Ofc 


x r — x s \)\ds(x r )ds(x s ) < 


II 

J J Oz, 


lQ k 

< 2ttR r R s 

< CR r Ra , 


In 2kR s \fr sin 


| d r — 0 S 


d0 r d0 s 


r0o 

/ ln(2fc R s y/r sin t)dt 

Jo 


where we have used integration by parts in obtaining the last inequality. Substituting 
the above estimate into (13TTTI we obtain 


f f ,, w \u s (x r ,X s )\ 2 , . , . 

// $(^,a; s )$(a; r , 2 )—r-- —ds(x r )ds(x s ) 

J jQ k ^ (£ r ,£ s ) 


<Ck- 2 (l + kd D ) 4 (kR s )-\ (3.12) 



























Next we estimate the integral in T r x r s \0/ e . Since t\H^ (t) | 2 is an increasing function 
of t > 0 |24l P- 446], we have for ( x r ,x s ) G T r x T s \Clk, \x r — x s \ > 1/(2 k), and thus 


\x r - X s \\u\x r ,x s )\ 2 > — k 1 


H t 


(i) 


= Ck~ 


which implies by using Lemma 13.21 and (|3.8|) again that 



xr s \n k 


$(z,x s )$(x r , 


z) 


\u ( x r , x s ) | 

U^Xr.Xs) 


ds(x r )ds(x s ) 


< Ck- 2 (l + kd D ) 4 (kRs)~ %/2 . 

(3.13) 


This completes the proof by combining the above estimate with (|3.12)) . □ 

Now we turn to the estimation of the third term at the right hand side of (USD. 
Denote S = ( kR s ) -1 / 2 and @<5 := {(0 r , 0 S ) G (0, 2tt) 2 : 1 0 r — 0 S ± m tt\ < S, m = 0,1, 2} 


and Qs •— {(^-r; ^s) £ T r x T s . ($ r ,$ s ) G @( 5 }. 
Lemma 3.5. We have 


! ff x s )Q(x r , z)u s (x r ,x s ) U .^ ri Xs ^ ds(x r )ds(x s ) 

J Jo 


u i {x r ,x s ) 


< C(l + kd D ) 2 (kR s y l ' 2 . 


Proof. The proof follows easily from Lemma 13.21 and (j3.8p and the fact that 

\Qs\ < CR r R s {kR s )-y 2 . □ 

To estimate the integral in T r x T s \Qs, we recall first the following useful mixed 
reciprocity relation [16], [23j P.40]. 

Lemma 3.6. Let = e 1 ^/VSirk. Then u^(d,x s ) = 7 m ^ s (x s ,—d) for any 
x s G M 2 \Z),d G S' 1 , where u^,(d,x s ) is the far field pattern in the direction d of the 
scattering solution of with u l (x ) = 4>(x,x s ) and u s (x,d) is the scattering 

solution of with the incident plane wave u l (x) = e lkx ' d . 

Lemma 3.7. Let u^(x s , —x r ) be the far field pattern in the direction x s of the 
scattering solution of the Helmholtz equation with the incident plane wave e ~ lkXr ' x . 
Then \u^(x 8 , -x r )\ + (kd D )~ 1 \du^(x s , -x r )/d0 8 \ < Ck~ x ! 2 (\ + kd D ) 2 . 

Proof The proof follows from the definition of the far field pattern in (|3.4|) with 
g(y) = —e~ lkXr ' y on Tp, Lemma l3Tl and ()3.1D . Here we omit the details. □ 

The following lemma is essentially proved in [ 8 ] Theorem 2.5]. 

Lemma 3.8. For any x G M 2 \D, the solution of the scattering problem (lT2l)-(lT3l) 
satisfies the asymptotic behavior: 


w(x) 


p ik\x\ 

—=w°°(x) + 'y(x), 
VFI 


where | 7 (x)| < C(l + kd D f(k\x\) 3/2 \\g\\m^(r D )- 

Proof First we have the following integral representation formula 


w(x) = / 
Jtt 


d$(x,y) dw(y) 

w(y) \ - ®{x,y)' 


du(y) 


dv{y) 


By the asymptotic formulae of Hankel functions 
1/2 


ds{y), Va; g R 2 \D. 
H P.197], for n= 1 , 2 , 


H n\t) = (Jf) + R n (t), \R n (t)\ < Ct~ 3 7 Vt > 0, 


(3.14) 


(3.15) 
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and the simple estimate \k\x — y | — k(\x\ — x • y)\ < C (k\y\) 2 (k\x\) 1 for any y E T^, 
x E r 2 \d, we have 

if p ife|x| 

= ^r=f y + 7o ^’ 

Vo7r/c y/\x\ 

d$(x,y) _ e ^ e ik\x\ g e -ikx-y 

dv{y) + ^ X ’ y >’ 

where \jo(x,y)\ + k~ 1 \ji(x, y)\ < C( 1 + k\y\) 2 (k\x\)~ 3 / 2 for some constant C inde¬ 
pendent of k and D. The proof completes by inserting ()3.16j) - ()3.17|) into ()3.14j) and 
using Lemma l3Tl and (|3.lj) . Here we omit the details. □ 

We also need the following slight generalization of Van der Corput lemma for the 
oscillatory integral |T2j P.152]. 

Lemma 3.9. For any —oo<a<b<oo, for every real-valued C 2 function u that 
satisfies \u'(t)\ > 1 for t E (a, b). Assume that a = xo < x\ < • • • < xn = b is a 
division of (a,b) such that u' is monotone in each interval (. Xi-\,Xi ) ; i = l,--- ,N. 
Then for any function <f defined on (a, b) with integrable derivative, and for any A > 0 ; 


(3.16) 

(3.17) 


r b 

/ e lXu ^(j)it)dt 

J a 


< (2N + 2)A" 


\4>(b)\ 


/ 


I <j>’{t)\dt 


Proof By integration by parts we have 

' g iA u(t) 


p 


e iX uW dt = 


iXu'(t) 


- [ e iXu ^ — 

Ja dt 


d ( 1 


iXu'(t) 


dt. 


Since u! is monotone in each interval (o^_i, xf), i = 1, • • • , N, and \u'(t)\ > 1 in (a, 6), 
we have 


f 

J a 


D iXu(t) _ 


d ( 1 


dt \iXu'(t) 


dt 


N 






f Xi A f 1 

J Xi _! dt \u'(t) 


dt 


< 2NX 


-l 


which implies | e lXu ^dt\ < (2 N + 2)A 1 . For the general case, we denote F(t) = 
ft e iA u(s) g s anc [ use integration by parts to obtain 

f <f(t)e lXu ^ dt = cj)(b)F(b) — f F(t)(/)' (t)dt. 

J a J a 

This completes the proof by using \F(t)\ < (2 N + 2)A -1 . □ 

Lemma 3.10. We have 


x s )$(x r , z)u s (x r , x s ) u Xs ) ds(x r )ds(x s ) 

u l [x r ,x s ) 


k 2 ff 

J JF r xF s \Qs 

< C( 1 + kd D ) 3 {kR s r 1 / 2 + C(1 + k\z\) 2 {kR s )-\ 


Proof We first observe that for ( x r ,x s ) E T r x T S \Q < 5 , we have 

Or -0 S 


k\x r — x s \ > 2 kR s \fr 


sm ■ 


5 1 

> 2fci? s \/rsin — > ~(kR s ) 1/2 y/r, 
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where we have used the fact that sin t > t/2 for t G (0, tv/ 2). Thus by ()3.15|) we obtain 

= e -™\*r-x .\+if +/Oo( x r ,x s ), (3.18) 

u [x r , X s ) 

where \po(x r ,x s )\ < C(kR s )~ 1 ^ 2 . Similar to (|3.16|) we have 

i— ikR 

Hz,x s ) = -^=J-j=e~ ik&s ' z + p 1 (z,x s ), (3.19) 

p if AkR r 

®(z,x r ) = ^/8^~^W e ~ ik&r ' Z + Pl ^ z,Xr ^ ( 3 - 20 ) 

where \pi(z,x s )\ < C(1 + k\z\) 2 (kR s )~ 3 / 2 , \pi(z,x r )\ < C{1 + k\z\) 2 (kR r )~ 3 ^ 2 . Next 
by Lemma 13.81 the mixed reciprocity in Lemma 13.61 and (|3.lj) we have 

e \kR r 

U s (x r ,x s ) = —=U™(x r ,X s )+p 2 (x r ,X s ) 

V Rr 

pi kR r 

= “7=7' mU S (x s ,-X r ) + p 2 (x r ,X s ) 

V R'p 

gik(R r -\-R s ) 

= - , -£ r ) + p 2 (x r ,x s ) + p 3 (x r ,x s ), (3.21) 

V R'p Rs 

where — x r ) is the far field pattern of the scattering solution of the Helmholtz 

equation with the incident plane wave u l = e -*kx r -x an q 

\p 2 (x r ,x s )\ < C( 1 + kd D ) 3 (kR r )~ 3/2 ||$(-,x s )|| ff i/ 2 (rD ) 

< (7(1 + kd D ) 4 {kR r )- 3/2 {kR s )~ 1/2 , 

\p 3 (x r ,x s )\ < (7(1 + fcd£)) 3 (fci? r .) _1/2 (A:i? s ) _3/2 ||e _l/c:i:r ' :i: || ff i/ 2 ( rD ) 

< (7(1 + kd D ) 4 (kR r )- 1 / 2 (kR s )- 3 / 2 . 

Combining (|3.18j) - (j3.21j) we have 


k 2 [[ $>(z, x s )&(x r , z)u s (x r ,x s ) U .[ Xr ' ds(x r )ds(x s ) 

JJr r xr s \Q 5 u\x r ,x a ) 


— k 2 R s R s II 3 >(z, x s )<f>(x r , z)u s (x r , x s )^r 


H . 

J J(0,27r) 2 \© (5 
_ 7m ^ 2\h(R r +Rs) 


U^Xr.Xs) 


d 0 r d 0 s 


87T 


(0,27t) 2 \© 5 


U^Xr.Xs) 

(j){ 6 r , 0 s )e~ 2 lk ^ Xr ~ Xs ^d 0 r d 0 s + p(z), ( 3 . 22 ) 


where (j)( 0 r ,Q s ) = —x r )e lk (xs+x r )-z an q py Lemma 13221 


|p(*)| < C (1 + kd D ) 2 (kR s )-^ 2 + C (1 + fcd D ) 2 (feii fl )- 1 / 2 + C (1 + a^|) 2 clr s )-\ 
< C (1 + yi)) 2 p s )- 1/2 + C{ 1 + fe| 2 f|) 2 (fei? a )^, 


where the second inequality follows from the fact that we are interested in the situ¬ 
ation when kR s is sufficiently large such that (1 + kd]j) 2 (kR s )~ 1 ' 2 <C 1 . Now direct 





















calculation shows that 


[[ 4>(0 r ,e s )e- 2ik ^- x ^de r d0 s 

J 0 , 27 r ) 2 \© 5 


[j 


(0 r +< s ),0 r .+7r—<5)U(0 r +7T+(5,0 r +27r — <5) 


0(<9 r ,6> s )e- 2i/c| ^-^!(i<9 s ^(9 r 


+ 

pTT — S 

[ 4>{9 r , 9 s )e~ 2lk ^ Xr ~ x ^d9 s d6 r 


J 8 J 

(0 ,6 r —<5)u(0 r -|-(5,0 r >-|-7r—<5)u(0 r +7r+(5,27r) 

+ 

P7T-\-S 

f 4>(0 r , 9 s )e~ 2ik ' [Xr ~ Xs ' [ d9 s d9 r 


Jtt-8 J 

( 0 ,. 7T • -V.o,. <5)U(« r t rf.fl,- • 7T <5} 

+ 

p2tt—8 

f (t>( y e ri e s )e- 2ik ^- x ^de s de, 


J 7T + <5 

J (0,0 r —7 T—<5)U(# r —7r+(5,0 r —<5)U(# r +<5,27r) 

+ 

r 

[ (f>{e r ,e s )e~ 2ik]Xr ~ xA de s de r 


to 

1 

(6 r —27T-\-8,0r—7r—8)\J(0 r —7T-\-8,6 r —8) 


: — Ii + • • * H- I5. 


(3.23) 


By Lemma 13771 and 5 = ( kR s ) -1 / 2 we have II 1 +I 3 +I 5 I < C/c _1 / 2 (l + fc(i^) 2 (/ci? s ) _1 / 2 . 

We will use Lemma 13.91 to estimate I 2 and I 4 . For that purpose, denote by 
v(6 s ) = —y/l + t 2 — 2t cos(# r — 0 S ). We have i/(0 s ) = r sin(0 s — 0 r )/v(6 s ) and thus 
\v f (0 s )\ > t\ sin S\/\v(6 s )\ > > 6/ 4 = ±(fciJ a )- 1/2 for (0 r ,0 a ) G T r x r fl \9«. 

Moreover, v"(0 s ) = — t 2 (cos(# s — 0 r ) — t)(cos(# s — 0 r ) — r _1 )/n(6 > s ) 3 which implies 
v'(6 s ) is piecewise monotone in (0, 27r) for any fixed 0 r G (0, 27t) since r > 1. 

Now since — 2\k\x r — x s \ = 2ifci2 s i;(0 S ), we obtain by Lemma [3791 and Lemma [3771 
that for 0 r G (S,7r — S), 


r0 r — (5 

Jo ' 


4>(0 r , 9 s )e~ 2ik ' [Xr ~ x ^d9 s 


< Ck~ 1/2 (1 + kd D ) 3 (kR s )~ 1/2 . 


The other integrals in I 2 and I 4 can be estimated similarly to obtain 

|I 2 | + M < Ck~^ 2 { 1 + kd D f{kR s )~ 1/2 . 

This completes the proof by (|3.22j) . □ 

The following theorem is the main result of this section. 

Theorem 3.1. For any z G Q, let 'ip(x^z) be the scattering solution to the 
problem: 

Ai/;(x, z) + k 2 ^p(x^ z) = 0 in M 2 \D, ^(x, z) = — Im 4>(x, z) on Fjj. (3.24) 


Then if the measured field \u(x r ,x s )\ = \u s (x r , x s ) + u l (x r , x s )| with u s (x,x s ) satis¬ 
fying the problem (HU-dLSl) with the incident field u l (x,x s ) = <L(x, x s ), we have 

I{z) = k f l^ 00 ^, z)\ 2 dx + Wj(z), \/z G fl, 

where \wj(z)\ < C( 1 + &dc>) 4 (&;ii ^) _1 / 2 + (7(1 + kdo) 2 {l + k\z\) 2 (kR s )~ 1 . 


Proof By (|3.7 p , Lemma 13.51 and Lemma 13.101 we are left to show that 


—k 2 Im 


4>(z, x s )4>(x r , z)u s (x r , x s )ds{x r )ds{xs) = k / |' 0 OO (^, 2 :)| 2 dx + £(z), 

./S ’ 1 
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with \({z)\ < (7(1 + kd]j) 2 ( 1 + kdjj + k\z\) 2 (kR s )~ 1 . This can be done by a similar 
argument as that in [3j Theorem 3.2]. For the sake of completeness, we include a 
sketch of the proof here. 

By the corollary of Helmholtz-Kirchhoff identity in |3j Lemma 3.2], for any z,y £ 

ft, 

k j <$>(x r , z)<$>(x r ,y)ds(x r ) = lm$(z,y) + w r (z,y), (3.25) 

where \w r (z,y)\ + k~ 1 \w r (z^y)\ < (7(1 + k\y\ + k\z\) 2 (kR r )~ 1 . Now by the integral 
representation formula 


U S (Xr 


,i-> ■ L ( 




<Mi/) 9u(y) 


have 


fr D 

where by m 


k J $(x r , z)u s (x r ,x s )ds(x r ) 

/ u s (y,x s ) 

>/r D L 


dlm $(z,y) du s (y,x s ) 


ICl(^,^s)l = 


dv{y) <My) 

-dw r (z,y) du s (y,x s ) 


Im ®(z,y) ds(y) + Ci(z,x 3 ), 


r ^ 

/ u s (y,x s )- 
J r n L 


/r D L dv{y) dv{y) 

< C(1 + Md) 2 (1 + Md + k\z\) 2 (kR s )~ 3 / 2 
By the definition of the imaging function I(z), we have then 

<9Im<l>(2:,y) dv s (y,z ) 


w r (z,y) ds(y) 


I(z) = — Im 


r 

/ v s (y,z)- 

JT d 1 


dv{y) dv{y) 


Im®(z,y) ds(y) +( 2 (z), (3.26) 


where v 8 (y,z) = k f r <&(x s , z)u s (y,x s )ds(x s ) and 

\( 2 {z)\ = k J §(x s ,z)Qi(z,x s )ds(x s ) 

< C{ 1 + kdn) 2 ^ + kdo + k\z\) 2 (kR s ) 
Taking the complex conjugate we get 


(3.27) 


3 (y,z) = kJ $(x s ,z)u s (y,x s )ds(x s ). 


Therefore, v s (y,z) can be viewed as the weighted superposition of v s (y. x s ). Then 
v s{y,z) satisfies the Helmholtz equation 

A v v 8 (y, z) + k 2 v s (y, z) = 0 in K 2 \D. 

On the boundary of the obstacle T^, we have 

v s (y, z) = k $(x s ,z)u s (y,x s )ds(x s ) 

Jr s 

= -k $(x s ,z)$(y,x s )ds(x s ) 

Jr„ 

=-Im $(z,y)-w 8 (z,y), Vy e T^, 
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where, similar to (13.251) , \w a {z,y)\ + k 1 \w s (z,y)\ < C(1 + k\y\ + k\z\) 2 (kR s ) 1 . By 
the definition of ip we have 


I(z) = -Im / ip(y,z) 
Jt d l 


<9Im $>(z,y) dip(y,z) 


dv{y) 


dv(y) 


Im $(z,y) ds(y) + ( 3 (z), 


= —Im 


L 


dip(y,z) 


'<P(y,z)ds(y) + C3(z). 


!Tn du(y) 

where we have used the boundary condition of ip on Td in the second inequality and 


&(z) = C, 2 {z) + Im 


/ <i>(y,z) 

Jr n L 


dlm $(z,y) d<j>(y,z) 


>r D L 9u(y) dv{y) 


■Im $>(z,y) ds(y). 


Here <p(y,z ) is the solution of scattering problem (|3.2|) - (|3.3|) with the boundary con¬ 
dition g = w s (z,y). Now by ((3.1 p and ()3.27| ) we then obtain 


|C3(^)| < C( 1 + kd D ) 2 (l + kd D + k\z\) 2 (kR s )~ 


Now the theorem is proved by using (1331) . □ 

We remark that ip(x, z) is the scattering solution of the Helmholtz equation with 
the incoming field J$(k\x — z\). It is well-known that Jo(t) peaks at t — 0 and decays 
like t -1 / 2 away from the origin. The source of the problem (13.24|) will peak at the 
boundary of the scatterer D and becomes small when z moves away from dD. Thus 
we expect that the imaging function I{z) will have a large contrast at the boundary 
of the scatterer D and decay outside the boundary dD. This is indeed observed in 
our numerical experiments. 

4. Extensions. In this section we consider briefly the imaging of the penetrable 
and impedance non-penetrable obstacles with phaseless data. We first consider the 
imaging of impedance non-penetrable obstacles with the phaseless data, in which 
case, the measured phaseless total field \u(x r ,x s )\ = | u s (x r ,x s ) + u l (x r ,x s ) |, where 
u s (x,x s ) is the radiation solution of the following problem: 


A u s + k 2 u s = 0 

in R 2 \D, 

(4.1) 

du s . . 

du l , \ ■ 


-1- ikr](x)u = 

av 

— —- ikrj(x)u on Tp. 

Ob' 

(4.2) 


Here rj(x) > 0 is the impedance function. The well-posedness of the problem (|4.1[) - 
((4.2 1) is well-known [7, [18]. By modifying the argument in section 3 and [3j Theorem 
3.2] we can show the following theorem whose proof is omitted. 

Theorem 4.1. For any z E Q, let ip(x, z) be the radiation solution of the problem 

Aip(x, z) + k 2 ip(x, z) = 0 in R 2 \D , 

+i= - 3Im ! (l ' z) -ih,wim*(*. 2 ) on r D . 

OV Ob' 

Then if the measured field \u(x r ,x s )\ = | u s (x r ,x s ) -\- u l (x r ,x s )\ with u s (x,x s ) satis¬ 
fying (|4.ip - (|4.2p . we have, for any z E Q, 

f(z) = k / \ip°°(x, z)\ 2 dx + k / rj(x) \ip(x, z) + Im$(r, z)\ 2 dx + Wf(z), 

Js 1 J r D 

ll 























where \wj(z)\ < C( 1 + kd]j) A (kR s )~ 1 / 2 + C( 1 + &dc>) 2 (l + /c|^|) 2 (fci? s ) _1 . 

For penetrable obstacles, the measured total field \u(x r ,x s )\ = \u s (x r ,x s ) + 
id(x r ,x s )|, where u s (x,x s ) is the radiation solution of the following problem 

A u s + k 2 n(x)u s = —k 2 (n(x) — 1 )u l (x,x s ) in M 2 (4.3) 

with n(x) G L°°(1R 2 ) being a positive function which is equal to 1 outside the scatterer 
D. The well-posedness of the problem under some condition on n(x) is known m- By 
modifying the argument in section 3 and in |3j Theorem 3.1], the following theorem 
can be proved. Here we omit the details. 

Theorem 4.2. For any z GO, let ip(x, z) be the radiation solution of the problem 

A'lp + k 2 n(x)^p = — k 2 (n(x) — l)Im 4>(x, z) in M 2 . (4.4) 

Then if the measured field \u(x r ,x s )\ = \u 3 (x r ,x s ) + u l (x r ,x s )\ with u s (x,x s ) satis¬ 
fying (|4.3p . we have 

I(z) = k f \fj°°(x, z)\ 2 dx + Wf(z) VzeQ, 

Js 1 

where \w f (z)\ < C{ 1 + kd D ) 4 (kR s )^ 2 + C{ 1 + kd D ) 2 (l + k\z\) 2 (kR,)- 1 . 

We remark that for the penetrable scatterers, ip(x,z) is again the scattering so¬ 
lution with the incoming field Im$(r,z). Therefore we again expect the imaging 
function I{z) will have contrast on the boundary of the scatterer and decay outside 
the scatterer. 

5. Numerical examples. In this section, we show several numerical experi¬ 
ments to illustrate the effectiveness of our RTM algorithm with phaseless data in 
this paper. To synthesize the scattering data we compute the solution u(x,x s ) of 
the scattering problem (HD-drai by standard Nystrom’s methods [8]. The boundary 
integral equations on T^ are solved on a uniform mesh over the boundary with ten 
points per probe wavelength. The boundaries of the obstacles used in our numerical 
experiments are parameterized as follows, where 0 G [0, 27r], 

Kite: x\ — cos (0) + 0.65 cos(2$) — 0.65, x 2 = 1.5 sin(0), 
p-leaf: r(0) = 1 + 0.2 cos(p0), 

Peanut: x\ = cos($) + 0.2cos(3$), X 2 = sin (6) + 0.2sin(3$), 
Rounded-square: x\ = cos 3 ($) + cos(0), X 2 = sin 3 (0) + sin(0). 

The sources x s , s = 1, • • • ,N S , and the receivers x r , x r = 1, • • • , 7V r , are uni¬ 
formly distributed on T s and T r , that is, x s = R s (cosO s , sin0 s ), 0 S = jf-(s — l),s = 
1, 2,..., N s , and x r = R r (cosO r , sin 6 r ),6 r = jf-(r — 1) + r = 1,2,..., 7V r , so that 
X r 7 ^ • 

Example 5.1. We consider the imaging of sound soft obstacles including a circle, 
a peanut, a kite and a rounded-square. The imaging domain is Q = (—3,3) x (—3,3) 
with the sampling mesh 201 x 201. The probe wave wavenumber k = Air, N s = N r = 
128, and R s = R, = 10. 

The imaging results are depicted in Figure l5Tl which show clearly that our imaging 
algorithm can find the shape and the location of the obstacles using phaseless data 
regardless of the shapes of the obstacles. 

Example 5.2. We consider the imaging of a 5-leaf obstacle with impedance 
condition g = 5, a partially coated obstacle with g = 5 in the upper boundary and 
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Fig. 5.1. Example 5.1: Imaging results by RTM imaging function with phaseless data. 

Top row: circle (left) and peanut (right); Bottom row: 5-leaf (left) and diamond (right). 


rj = 1 in the lower boundary, a sound hard, and a penetrable obstacle with n(pc) = 0.25. 
The imaging domain is Q = (—3,3) x (—3,3) with the sampling grid 201 x 201. The 
probe wave wavenumber k = An, N s = N r = 128, and R s = R r = 10. 

Figure [521 shows the imaging results which demonstrate clearly that our imaging 
algorithm works for different types of obstacles without using any a prior information 
of the physical properties of the obstacles. 

Example 5.3. We consider the stability of the imaging function with respect 
to the additive Gaussian random noises using the phaseless data. We introduce the 
additive Gaussian noise as follows (see e.g. W: 

|^| noise = |ll| T ^noise: 


where \u\ is the synthesized phaseless total field and v no i S e is the Gaussian noise with 
mean zero and standard deviation pi times the maximum of the data \u\, i.e. z/ no i se = 
/imaxims, and £ ~ A/"(0,1). 

For the fixed probe wavenumber k = An, we choose one kite and one 5-leaf in 
our test. The search domain is Q = (—5,5) x (—2,4) with a sampling 201 x 201 
mesh. We set R s = 10, R r = 20, and N s = N r = 256. Figure 15.31 shows the 
imaging results for the noise level pi = 10%, 20%, 30%, 40% in the single frequency 
data, respectively. The imaging results can be improved by superposing the multi¬ 
frequency imaging result as shown in Figure 15.41 The left table in Table 15.11 shows 
the noise level, where a = £<max Xr)Xs |u(a; s ,av)|, \\u\\% = \u(x s .x r )\ 2 , 

■ T,s,r=l \ V 


2 - - -J—\" Ns ’ Nr \u ise (x s ,x r )| 2 . 


'noise 11 ,£2 — 


N s N r 
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Fig. 5.2. Example 5.2: In the top row, a sound hard 5-leaf obstacle (left) and a non-penetrable 
obstacle with the impedance rj = 5 (right). In the bottom row, a partially coated obstacle with 
= 5 on the upper boundary and rj = 1 on the lower boundary (left) and a penetrable obstacle with 
(x) = 1/4 (right). 



Fig. 5.3. Example 5.3: The imaging results using single frequency data added with additive 
Gaussian noise ji = 10%, 20%, 30%, 40% from (a) to (d), respectively. The probe wavelength is 
A = 0.5 and the sampling number is N s = N r = 256. 
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(c) 


(d) 


Fig. 5.4. Example 5.3: The imaging results using multi-frequency data added with additive 
Gaussian noise /a = 10%, 20%, 30%, 40% from (a) to ( d), respectively. The probe wavelengths are 
given by X = 1/1.8,1/1.9,1/2.0,1/2.1,1/2.2 and the sampling number is N s = N r = 256. 


R 

a 

Kll < 2 

11 ^noise 11 £ 2 

0.1 

0.002859 

0.013054 

0.002863 

0.2 

0.005717 

0.013054 

0.005708 

0.3 

0.008576 

0.013054 

0.008572 

0.4 

0.011435 

0.013054 

0.011424 


R 

a 

\\ u \\p 

11 ^noise 11 £ 1 2 3 4 5 6 7 8 9 

0.1 

0.003004 

0.013017 

0.003007 

0.2 

0.006009 

0.013017 

0.005996 

0.3 

0.009013 

0.013017 

0.008964 

0.4 

0.012018 

0.013017 

0.012008 


Table 5.1 

Example 5.3: The noise level in the case of single frequency data (left) and multi-frequency 
data (right). 
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